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In this review, we discuss the decoherence and thermalization of a quantum spin system interact- 
ing with a spin bath environment, by numerically solving the time-dependent Schrodinger equation 
of the whole system. The effects of the topologic structure and the initial state of the environment on 
the decoherence of the two-spin and many-spin system are discussed. The role of different spin-spin 
coupling is considered. We show under which conditions the environment drives the reduced density 
matrix of the system to a fully decoherent state, and how the diagonal elements of the reduced 
density matrix approach those expected for the system in the microcanonical or canonical ensemble, 
depending on the character of the additional integrals of motion. Our demonstration does not rely 
on time-averaging of observables nor does it assume that the coupling between system and bath is 
weak. Our findings show that the microcanonical distribution (in each eigenenergy subspace) and 
canonical ensemble are states that may result from pure quantum dynamics, suggesting that quan- 
tum mechanics may be regarded as the foundation of quantum statistical mechanics. Furthermore, 
our numerical results show that a fully decoherent quantum system prefers to stay in an equilibrium 
state with a maximum entropy, indicating the validity of the second law of thermodynamics in the 
decoherence process. 

PACS numbers: 05.30.-d 03.65.-w 03.65.Yz 05.45.Pq 

Keywords: Quantum Statistical Mechanics, Canonical Ensemble, Time-dependent Schrodinger Equation, 

Decoherence, Second Law of Thermodynamics 



Contents 

I. Introduction 
II. General Theory 

III. Quantum Spin System and Numerical 
Method 

IV. Decoherence of a Two-Spin System 

A. A Two-Spin System 

B. Decoherence Without Energy Dissipation 

C. Decoherence With Energy Dissipation 

D. Summary 

V. Evolution to the Ground State of a 
Two-spin System 

A. Symmetrical Coupling 

B. Antisymmetric Coupling 

C. Summary 

VI. Decoherence of a Many- Spin System 

A. Origin of the Microcanonical Distribution 

B. Summary 

VII. Thermalization of a Many-Spin System 

A. Origin of the Canonical Ensemble 

B. Summary 

VIII. Conclusion and Discussion 
IX. Acknowledgement 
References 



1 



a 



m 



m 

m 

12 

a 

i 

m 

m 



I. INTRODUCTION 

The manner in which a quantum system becomes effec- 
tively classical is of great importance for the foundations 
of quantum physics. It has become increasingly clear that 
the symptoms of classicality of quantum systems can be 
induced by their environments^. 

Intuitively, we expect that by turning on the interac- 
tion between the quantum system and the environment, 
the fluctuations in the environment will lead to a reduc- 
tion of the coherence in the quantum system. This pro- 
cess is called decoherence^r— . The existence of decoher- 
ence in the quantum system represents a challenge for 
the realization of quantum computation and quantum 
information processing^—, which are expected to rely 
heavily on quantum coherence. In general, there are two 
different mechanisms that contribute to decoherence. If 
the environment is dissipative (or coupled to a dissipa- 
tive system), the total energy is not conserved and the 
whole system relaxes to a stationary equilibrium state, 
for instance the thermal equilibrium state. In this review, 
we exclude this class of dissipative processes and restrict 
ourselves to closed systems in which a small quantum 
system is brought in contact with a much larger environ- 
ment. Then, decoherence is solely due to the fact that the 
initial product state (wave function of the quantum sys- 
tem times wave function of the environment) evolves into 
an entangled state of the whole system. The interaction 
with the environment causes the initial pure state of the 
quantum system to evolve into a mixed state, described 
by a reduced density matrix 1 ^, obtained by tracing out 
all the degrees of freedom of the environment ^ 11 : 17 : 18 . 

The decoherence programme is supposed to explain the 



macroscopic quantum superposition ( "Schrodinger cat" ) 
paradox, that is, the inapplicability of the superposition 
principle to the macroworld. The states that are "ro- 
bust" with respect to the interaction with the environ- 
ment are called pointer states 9 . If the Hamiltonian of 
the quantum system H$ is a perturbation, relative to 
the interaction Hamiltonian Hi nt , the pointer states are 
eigenstates of Hi n ®~ 11 . In this case, the pointer states 
are essentially "classical states" , such as states with def- 
inite particle positions or with definite spin directions 
of individual particles for magnetic systems. In general, 
these classical "Schrodinger cat states" , being a product 
state of individual particles forming the system, are not 
entangled. 

On the contrary, if the interaction Hamiltonian Hi nt is 
a perturbation, relative to the Hamiltonian of the quan- 
tum system H$, the pointer states are eigenstates of 
-ffs&ii. In this case, the pointer states are not necessary 
classical-type states, they may be "quantum" states such 
as standing waves, stationary electron states in atoms, 
tunncling-split states for a particle distributed between 
several potential wells, singlet or triplet states for mag- 
netic systems, etc—. This may explain, for example, 
that one can observe linear atomic spectra - the initial 
states of an atom under the equilibrium conditions are 
eigenstates of its Hamiltonian and not arbitrary superpo- 
sition thereof. Some less trivial pointer states have been 
found in computer simulations of quantum spin systems 
for some range of the model parameter o 20 i 21 . In fact, the 
evolution to equilibrium of quantum spin systems is still 
an open issue. Recently, the effect of an environment 
of JV > 1 spins on the entanglement of quantum spin 
systems has attracted much attention 2 ^—. 

Furthermore, recent simulation results of quantum 
spin systems show that different statistical ensembles 
such as the microcanonical (per eigenenergy subspace) 
and the canonical ensemble, could arise from the distri- 
bution of these pointer states (eigenstates) in the mixed 
stated. This direct connection of quantum dynamics and 
statistical mechanics gives an explanation of a basic pos- 
tulate in statistical mechanics: a generic "system" that 
interacts with a generic environment evolves into a state 
described by the canonical ensemble. Experience shows 
that this is true but a detailed understanding of this pro- 
cess, which is crucial for a rigorous justification of sta- 
tistical physics and thermodynamics, is still lacking^ 3 - - — . 
In particular, in this context the meaning of "generic" is 
not clear. The key question is to what extent the evo- 
lution to the equilibrium state depends on the details of 
the dynamics of the whole system. 

Earlier demonstrations of the fact that the system can 
be in the canonical ensemble state are based on the Er- 
godic averages, that the expectation values of the dynam- 
ical variables of the system approach their values for the 
subsystem that is in the thermal equilibrium state^r— , 
or do not consider the dynamics of the system but assume 
that the state of the whole system has a special property 
called "canonical typicality"— ~—. There are two basic 



assumptions in the derivation of the canonical typical- 
ity, one is that the whole system is in the microcanoni- 
cal ensemble, another is that the interaction between the 
system and the environment is so small that it can be 
neglected. The theory of the canonical typicality is ki- 
netic rather than dynamic, and it is yet unclear under 
which conditions the whole system will evolve into the 
region in Hilbert space where its subsystems are in the 
thermal equilibrium state— A very different setting to 
study nonequilibrium quantum dynamics is to start from 
an eigenstate of some initial Hamiltonian and push the 
system out of this state by a sudden change of the model 
parameters^—. To the best of our knowledge, it has not 
yet been shown that this approach leads to the establish- 
ment of the canonical equilibrium distribution. We also 
want to draw attention to the fact that a demonstration 
of relaxation to the canonical distribution requires a sys- 
tem with at least three different eigenenergies because a 
diagonal density matrix of a two-level system can always 
be represented as a canonical distributio n 61 ^ 2 . 

In this review, we focus on recent results concern- 
ing the decoherence and thermalization in quantum spin 
system s 20 ' 21 ' 23 " — . In section II, we give a general the- 
ory about decoherence and thermalization, and intro- 
duce several quantities that measure the effect of deco- 
herence and thermalization. In section III, we introduce 
the quantum spin model and the methods used in the 
numerical simulation. In section IV, we focus on the de- 
coherence of a two-spin quantum system, and consider 
the effect of different topological structures and different 
types of spin-spin coupling. In section V, we discuss how 
a Hcisenberg two-spin system evolves to the ground state. 
In section VI, we consider the decoherence of many-spin 
systems. We show under which conditions the environ- 
ment drives the reduced density matrix of the system to 
a fully decoherent state, which is described by the mi- 
crocanonical distribution per eigenenergy subspace. In 
section VII, by introducing the energy dissipation of the 
many-spin system, we show that the diagonal elements 
of the reduced density matrix approach those expected 
for the system in the canonical ensemble. Section VIII 
contains the conclusion and a brief discussion about the 
second law of thermodynamics in quantum systems. 



II. GENERAL THEORY 

In general, the state of a closed quantum system is 
described by a density matrb i 65 i 66 p. The state of a 
quantum system interacting with an environment is rep- 
resented by the reduced density matrix p{t), obtained 
by tracing out all the degrees of freedom of the environ- 
ment. Decoherence of the quantum system means that 
the amplitude of the off-diagonal terms in the reduced 
density matrix become smaller, and full decoherence cor- 
responds to all off-diagonal terms being zero. Here we 
consider the case that the interaction Hamiltonian Hi nt 
is a perturbation, but not necessary very small, relative 



to the Hamiltonian of the quantum system Hs, and the 
reduced density matrix of the quantum system is repre- 
sented in its energy eigenstates. The decoherence of the 
quantum system can also be monitored by the time de- 
pendence of the quadratic entropy Ss (t) = 1 — Trp 2 (£) 
and the Loschmidt echo^ L (t) = Tr (p (t) po (£)), where 
po (t) is the density matrix for Hi nt = 0. 

The microcanonical ensemble is a mixed state where 
all accessible eigenstates have equal weight. The micro- 
canonical distribution per eigenenergy subspace is char- 
acterized by a density matrix that is diagonal with re- 
spect to the eigenstates of the system Hamiltonian, and 
the diagonal elements which belong to the degenerate en- 
ergy eigenstates are equal. A state with microcanonical 
distribution per eigenenergy subspace is a microcanonical 
ensemble if it has only one accessible eigenenergy. 

The canonical ensemble is a mixed state where the di- 
agonal elements take the form exp(—{3Ei), (3 = 1/ksT is 
proportional to the inverse of the temperature T (JtB is 
Boltzmann's constant) and the E^s denote the eigenen- 
ergies. 

The distribution of the state of a quantum system is 
the microcanonical or canonical ensemble only if it is in 
a fully decoherent state. 

The time evolution of a closed quantum system is 
governed by the time-dependent Schrodinger equation 
(TDSE)2£i££. If the initial density matrix of an isolated 
quantum system is non-diagonal, then, according to the 
TDSE, its density matrix remains nondiagonal and never 
approaches the thermal equilibrium state with the canon- 
ical or microcanonical distribution. Therefore, in order to 
thermalize the system S, it is necessary to have the sys- 
tem S interact with an environment (E), also called the 
heat bath. Thus, the Hamiltonian of the whole system 
(S + E) takes the form H = H S +H E + H SE , where H s 
and He are the system and environment Hamiltonian, 
respectively, and Hse describes the interaction between 
the system and the environment. 

The state of the system S is described by the reduced 
density matrix 



p(t) = Tr eps+e it) 



(1) 



where ps+E (t) is the density matrix of the whole system 
at time t and Tr^; denotes the trace over the degrees of 
freedom of the environment. 

The coherence of the system is conveniently character- 
ized by a(t), which is a global measure for the size of the 
off-diagonal terms of the reduced density matrix, defined 
by 



a(t) 



\ 



E E M*)i s 



(2) 



i—l j—i-\-l 



Here n denotes the dimension of the Hilbert space of 
system S and py(t) is the matrix element (i,j) of the 
reduced density matrix p in the representation that di- 
agonalizes Hs- If <r{t) =0 the system is in a state of full 
decoherence. 



The difference between the state p (t) and a micro- 
canonical distribution per eigenenergy subspace can be 
characterized by a(t) (measure of decoherence) and 7(f) 
(measure of the difference between the diagonal terms 
corresponding to the degenerate eigenstates), which is 
defined by 



lit) 



\ E E M*) " PoM 2 $ (Ei - E 3 ). (3) 
\ i=i j=i+i 



The system is in the microcanonical distribution per en- 
ergy subspace if and only if a{t) — and j(t) = 0. 

The system S is in its thermal equilibrium state only 
if the reduced density matrix takes the form 



e-^/Tr 



se 



-PH S 



(4) 



where Tr,g denotes the trace over the degrees of freedom 
of the system S. The difference between the state p (t) 
and the canonical distribution p is represented by a{t) 
and 5(t), defined by 



(5) 
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ith 
hit) 


n / n \ 2 

E U*w- e ~ b(t)B 7E e_Kt)B * > 

i=l V i=l ) 
<j,E^E 3 i^Pnit) ~ fopii(t)]/(Ej - E t ) 


0(1) 


Y.iKj^i^E, l 



(6) 



As the system relaxes to its canonical distribution both 
a(t) and S(t) vanish, b(i) converging to /3. 



III. QUANTUM SPIN SYSTEM AND 
NUMERICAL METHOD 



Most theoretical investigations of decoherence have 
been carried out for oscillator models of the environment 
for which powerful path-integral techniques can be used 
to treat the environment analytically 1 ' ?, ' ls . On the other 
hand, it has been pointed out that a magnetic environ- 
ment, described by quantum spins, is essentially different 
from the oscillator model in many aspects^. For the sim- 
plest model of a single spin in an external magnetic field, 
some analytical results are known^. For the generic case 
of two and more spins, numerical simulation— _ — is the 
main source of theoretical information. 

A generic quantum spin model can be described by the 
Hamiltonian H — Hs + He + Hse where 



H s = - 



H f 



H : 



SE 



Ns-1 N s 

E E E J tiS?sf, (7) 

i— 1 j—i-\-l a—x.y,z 

N-l N 

E E E ^^< (8) 

i— 1 j—i-\-l a=x,y,z 

N s N 

EE E A&s?jy. o) 

i— 1 j — 1 a— x,y,z 




FIG. 1: A typical configuration of a quantum spin system 
surrounded by a quantum spin bath (Ns — 4 and N = 18). 
The quantum spin system consists of four spin 1/2 particles, 
with orientation spin-up or spin-down. The bath spins are 
in a complicated random superposition state for which the 
expectation values (S x ) , {S y ) and {S z ) of each spin are all 
zero. 



Here the S a 's and I a, s denote the spin-1/2 operators of 
the system and environment respectively (we use units 
such that h and ks are one). An analytic solution of the 
TDSE can only be obtained for very special choices of the 
exchange integrals J"j, fi"- and Afj but it is straight- 
forward to solve the TDSE numerically for any choice of 
the model parameters. 

The state, that is the density matrix p(t) of the whole 
system at time i, is completely determined by the choice 
of the initial state of the whole system and the numerical 
solution of the TDSE. In this review, the initial state 
of the whole system (S+E) is a pure state. This state 
evolves in time according to 

ns n 

|*(t))=e^|*(0))=£5>(t,p,«)M, 

i— 1 p— 1 

where the states {| «,£>)} denote a complete set of or- 
thonormal states. In terms of the expansion coefficients 
c(i,p, t), the reduced density matrix reads 

n n 

P(t)i,j = Tr E ^2^2c*(i,q,t)c(j,p,t)\j,p)(i,q\ 

J>=1 q=l 



^2c*(i,p,t)c(j,p,t), 
P =i 



(10) 



which is easy to compute from the solution of the TDSE. 
We monitor the effects of decoherence by computing the 



matrix elements of the reduced density matrix p(t) of the 
quantum system. As explained earlier, in the regime of 
interest |A| <C \J\, the pointer states are expected to be 
the eigenstates of the quantum systems. Hence, we com- 
pute the matrix elements of the density matrix in the 
basis of eigenvectors of the quantum system. We also 
compute the time dependence of the quadratic entropy 
S s {t) = 1 - Trp 2 (t) and the Loschmidt echo^L (t) = 
Tr (po (t) p (t)) , where po(t) is the density matrix for 
Hse = 0. Another quantity of interest that can be ex- 
tracted from the solution of the TDSE is the local density 
of states (LDOS) 



LDOS(£0 



1 

2^ 



D 



+oo 



dt e-** t <tf(0)|e- li "|tt(0)) 



£|<*(0)|¥>*>| 2 *(25-25*), 



(11) 



fc=i 



where D = n + ns, {\(pk)}, and {E^} denote the dimen- 
sion of the Hilbert space, the eigenstates and eigenvalues 
of the whole system, respectively. The LDOS is "local" 
with respect to the initial state: It provides information 
about the overlap of the initial state and the eigenstates 
oi H. 

The notation to specify the initial state is as follows: 

(1) \GROUND) s is the ground state or a random su- 
perposition of all degenerated ground states of the sys- 
tem; 

(2) \GROUND) g is the state which has the energy 
close but not equal to that of the ground state; 

(3) \UU) s is a state in which all spins of the system 
are up, meaning that in this state, the expectation value 
of each spin is one; 

(4) \UD) S is a state in which two nearest- neighbor- 
spins of the system are antiparallel implying that in this 
state, the correlation of their z-components is minus one; 

(5) \UD) S is a state close to \UD) S , but the correlation 
of their z-components is larger than minus one; 

(6) \RR) S denotes the product state of random su- 
perpositions of the states of the individual spins of the 
system; 

(7) [RANDOM) o denotes a random superposition of 
all possible basis states. 

The same notation is used for the spins in the environ- 
ment, the subscript S being replaced by E. 

As we report results for many different types of spin 
systems it is useful to introduce a simple terminology 
to classify them according to symmetry and connectiv- 
ity. The terms "XY", "Hciscnberg" , "Heisenberg-type" , 
"Ising" , "Ising-type" and Ising± system refer to the fol- 
lowing cases: 

(1) XY: J*. 

(2) Hcisenberg: J? ■ 



- Jfj = J and J?j = 0; 
- T v - I z - T- 



(3) Heisenberg-type: Ji t] 
range [- | J\ , |J|]; 



are uniform random in the 



(4) Ising: J? 



jy 



and Jf 



J: 



Jf 



J v 



and Jf a are uniform 



*>j 



are random — \J\ 



(5) Ising-type: „ tJ 
random in the range [— \J\ , | J|]; 

(6) Ising± : Jf^ — Jfj = and Jf 
ox\J\. 

The same terminology is used for the Hamiltonian He 
of the environment and for the interaction Hamiltonian 
Hse- In our model, all the spins of the system interact 
with each spin of the environment. To characterize the 
connectivity K of the spins within the system (environ- 
ment), we use the term "ring" (K = 2) for spins forming 
a one-dimensional chain with nearest-neighbor interac- 
tions and periodic boundary conditions, "square-lattice" 
(K = 4) or "triangular-lattice" (K = 6) if the spins are 
located on a two-dimensional square or triangular lat- 
tice with nearest-neighbor interactions, and "spin glass" 
(K = N — 1) when all the spins within the system (envi- 
ronment) interact with each other with a Heisenberg-type 
interaction 68,69 . 

The time evolution of the whole system is obtained 
by solving the TDSE for the many-body wave function 
\ty(t)), describing the system plus the environment^r— . 
The numerical method that we use is described in Ref>22. 
It is based on the numerically exact Chebyshev poly- 
nomial decomposition of the operator U (t) = e~ ltH . It „£ _ 



is very efficient and conserves the energy of the whole 
system to machine precision. It is widely used in the 
simulation of quantum spin systems, such as the study 
of decoherenc o 20 ' 21 i 23 ~ — , the modeling of a quantum 
computer 8 -, the propagation of the quantum spin wave 86 , 
and the study of stability of the quantum domain wall2£, 
etc.. 

The simulation of the wave function is performed in 
the spin-up and spin-down basis, and it is convenient to 
use this basis to apply the Hamiltonian operator on the 
wave vector, which is a basic operation in the Chebyshev 
polynomial algorithm. Based on the wave function of the 
whole system in the up-down basis, we first calculate the 
reduced density matrix of the quantum system by tracing 
out the degrees of freedom of the environment, then di- 
agonalize the Hamiltonian Hs to get all the energy eigen- 
values and the corresponding eigenstates of the quantum 
system, and finally transform the reduced density matrix 
from the up-down basis to the energy eigenstate basis. 
With the reduced density matrix we can calculate all the 
physical quantities of the quantum system, such as the 
energy, the entropy, the measure of decoherence <r, the 
measure of the distribution 7, 5 and the effective inverse 
temperature b, etc. 



IV. DECOHERENCE OF A TWO-SPIN SYSTEM 

A. A Two-Spin System 

The quantum state of a two-spin system is completely 
determined by its reduced 4x4 density matrix. Although 
the reduced density matrix contains all the information 
about the quantum system, it is often convenient to char- 



acterize the state of the two-spin system by other quan- 
tities such as the correlation functions (Si • S 2 ), (S1S2), 
and (SfS'2), the single-spin magnetizations (Sf), (5*1), 
and M = (Sf + Sf), and the concurrence C 8 ^. In Ta- 
ble HI we show the values of these quantities for different 
states of the two-spin system. 

The concurrence, which is a convenient measure for the 
entanglement of the two spins is defined a a 89 ' 90 



C (p) = max(0, Ai - A2 - A3 - A 4 ), 



(12) 



where the Aj are the eigenvalues, in decreasing order, of 
the Hermitian matrix 



R - vVppVp- 



(13) 



Here p is the reduced density matrix of central spin pairs 
based on the standard basis |t) It), It) II), 1 1) It), It) It), 
and 



where a„ 



p = (a y ® a y )p*(a v <g> a y ), 



-i 



(14) 



„ I and p* is the complex conjugate 



In fact, the concurrence C is a measure between the 
state \ip) and the state with the two spins flipped i[) 



c = |(v#) 



(15) 



The singlet state, \ip) = (|f) |t) - It) It)) /"v/2 is un- 
changed under flipping two spins, therefore C = 1. 
The triplet state \ifj) = (|f) |t) + It) It)) /V% is also un- 
changed under flipping two spins, so C = 1. For |t) It), 
It) It), It) It), and |t) It), the state is totally different if 
the two spins flip and then C = 0. 

In the case that the two spins are coupled by the 
isotropic Heisenberg interaction, the Hamiltonian of the 
system simplifies 

H s = -JSiSa, (16) 

and the four eigenstates of Hs are given by 

\Ti) = ltt) = |l), 
\S)= ^-^=12), 



V2 
Itt) + Itt) 

|T_0 = HI) = |4) , 



\T n ) = 



= 13), 



(17) 



satisfying 

H s \S) = E s \S) , H s |Ti, ,-i) = E T |T 1A -i) , (18) 

where E s = 3 J/4 and E T = -J/4. 

From Table HI it is clear that the singlet state \S) is the 
most easily distinguished state as the two-spin system is 



TABLE I: The values of the correlation functions (Si • S2}, {Sf Sf), {S1S2), the total magnetization M, the concurrence C and 
the magnetization (Si) for different states of the two-spin quantum system. 



\<p) 


(Si ■ S 2 ) 


(SfSS) 


(S?Sf> 


M 


c 


(S!) 


(\n)-m))/V2 


-3/4 


-1/4 


-1/4 





1 





(\n) + \w)/V2 


1/4 


-1/4 


1/4 





1 





(\W-\U))/V2 


1/4 


1/4 


-1/4 





1 





(Itt) + |U» A/2 


1/4 


1/4 


1/4 





1 





in) 


-1/4 


-1/4 











1/2 


i;t> 


-1/4 


-1/4 











-1/2 


in) 


1/4 


1/4 





1 





1/2 


iu) 


1/4 


1/4 





-1 





-1/2 



in the singlet state if and only if (Si • S2) = —3/4. To 
identify the other states, we usually need to know at least 
two of the quantities listed in Table [I] For example, to 
make sure that the system is in the triplet state \Tq), 
the values of (Si ■ S2) and (Sf 5|) should match with the 
corresponding entries of Table HI Likewise, the two-spin 
system will be in the state \\f) if (Si • S 2 ) = 1/4 and 
M = 1. 

If the interaction between the Heisenberg quantum sys- 

tern and the environment is isotropic, that is, if A^ - J = 
r all i,j, then the Hamiltonian Hse 



A (y) = A {z) = A fo 
is simplified as 



N 



Hsi 



■A(Si 



5> 



(19) 



which leads to [Hs,Hse] = 0. As shown in Refill, if 
the energy of the quantum system is conserved, then the 
decoherence process is determined by Hse, He, the ini- 
tial state of whole system |\&(io)), and the eigenstates 
of the quantum system. In other words, in this case, 
L (t) and \p(t)\ do not depend on J, which means that 
the relative value of A/ J has no effect on the decoher- 
ence process. Furthermore, if we take the interactions 
between the environment spins to be isotropic, that is, 
$7jj for all i,j, then the Hamilto- 



,(*) 



•>(y) 



nian 



Q 



» 



He 



AT-l N 

E E 



* *i,j*-i 



(20) 



commutes with Hse , and therefore He has also no effect 
on the decoherence process. 

In fact, since [Hs,Hse] = 0, the time evolution op- 
erator of the whole system e~ tHt can be represented 
as e -iHst e -i{H SE +H E )t_ The initial state of the quan _ 

turn system can be represented as |y(io)) = ^2k a k \fy, 
where {|A;)} and {Ek} are the eigenstates and corre- 
sponding eigenvalues of the quantum system, that is, 
H s \k) = E k \k). 

For an isolated quantum system (Hse = 0), the time 
evolution of the density matrix of the quantum system is 



given by 



Mi)=]Te^ (Bfc - £!)t a fc ar|fc)(Z|. (21) 



k.l 



If the quantum system is coupled to a bath with initial 
state (j)(to), the state of the whole system at time t is 
given by 



|*(t)> - e'^ |*(*o)> 

?- iE kt n , p -i(H SE +H E )t 
k 



^e- tJikt a k e~ l{HsE+ME>t \k) \cf>(t )) .(22) 



As [H S ,H SE } = 0, we have H SE \k) \<p(t )) 
\k) Mk \<j) (to)) and hence, the state at time t becomes 



|*(t))=^a fc e-^ t |fc)|^(t)) 



(23) 



where 

\4> k (t)) = e-^ M " +H ^ t \cf > (t )). (24) 

The density matrix ps+E (t) of the whole system is 

p s+E (t) - mt)) <*(*)i 



E- 

fej 



-i(E k -Ei)t n * 



a k at\k)\(t>k(t))(l\(<t>i(t)\ 



(25) 



and the reduced density matrix p (t) of the quantum sys- 
tem is 

p(t) = Tr E ps+E(t) 



E« 

k,l 



-i(E k -E l )t ni * 



a k at {4>i (t) \<t> k (*)) \k) {l\ 



(26) 



The Loschmidt echo L (t) of the quantum system can be 



calculated as 



L(t) = Tr(p(t)p (t)) 

k,l 



Tr^e-^-^aua; (ft (t) \cj> k (t)) \k) (l\ 



xJ2e- KE -- E - )t a m a* n \m)(n\ 

J2\a k \ 2 \ ai \ 2 (fr(t)\^(t)). 
k,i 



(27) 



It is clear that if [Hs, Hse] = 0, the decoherence pro- 
cess is determined by the initial state of the quantum 
system {ak} and the time evolution of {\<j>k (£))}• As 
shown in Eq. (j2~4")l . the {\4>k (t))} are determined by the 
initial state of the bath (f>(to), the eigenstates {\k)} of the 
quantum system, and the Hamiltonians Hse and He- 
The eigenvalues {E^} have no effect on the decoherence 
process. Thus, multiplying Hs by a constant does not 
change L (t) and the diagonal elements of the reduced 
density matrix ps (t). The time evolution of the absolute 
value of the off-diagonal elements 



lw(*)ttl = M*lvM*) !&(*)) 



(28) 



is independent of Hs- This means that the relevant val- 
ues of the coupling between the system spins (J) and 
between the system and the environment (A) have no 
effect on the decoherence process. 

If [H S ,H SE } = and [H SE ,H E ] = 0, then, Eq. flU 
becomes 



!&(*)> = e- iM ^e- iH ^\4> (to)), 



(29) 



and therefore we have 

(&(*)!&(*)> = (<Mio)K' i(Mfc - Ml)t l<Mio)>, (so) 

implying that \pg (t) M \ and L (t) do not depend on He- 
Our goal is to find under which conditions the quantum 
system can evolve into a classical mixed state, that is, the 
elements in the reduced density matrix satisfying 



\ P s(t) a \ = Oiffc^Z, 

\ps(t) H \ ? onk = i. 



(31) 



B. Decoherence Without Energy Dissipation 

In previous work22r— , it was shown that a frustrated 
environment, such as described by a Heisenberg-type He, 
can enhance the decoherence of the quantum system. 
The typical results of the (full) decoherence in the Heisen- 
berg two-spin system without energy dissipation are the 
following. 

In Fig. [5] and Fig. [3J we show the time evolution of the 
elements of the reduced density matrix p (t) for different 
connectivity K with same f2, or different Q with same K, 
for the case [Hs, Hse] = 0. 
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FIG. 2: Time evolution of the real part of the off-diagonal el- 
ement P23 (left panel) and the diagonal elements pn, . . . , pa 
(right panel) of the reduced density matrix of a Heisenberg 
two-spin system (J = —5), coupled via an isotropic Heisen- 
berg interaction Hse (A = —0.075 ) to a Heisenberg-type 
environment He (fi = 0.1) with different connectivity: (a) 
K = 0; (b) K = 2; (c) K = 4; (d) K = 6; (e) K = N - 1. 
Adapted from Ref.— . 



If | A | ^> VL\/K, in the absence of interactions between 
the environment spins (Q,yK = 0) and after the initial 
decay, the quantum system exhibits long-time oscillations 
(see Fig. HJa) (left)). As shown in Ref322, in the limit 
of a large environment (JV — > oo) 



Re P23 (t) 



1-bt 2 



cosujt, 



(32) 



where b = 7VA 2 /4, c = 6/2 and w = J- A. Equation (|3"2"|) 
clearly shows the two-step process, that is, after the ini- 
tial Gaussian decay of the amplitude of the oscillations, 
the oscillations revive and their amplitude levels of by a 
factor of 1/3 (see Refj^-). Due to conservation laws, this 
behavior does not change if we introduce an isotropic 
Heisenberg Hamiltonian in the environment (Q^ • = il 
for all a, i and j), independent of K . This is also con- 
firmed by our numerical results (not shown). 

If |A| w Q\/K~, the initial Gaussian decay of the 



Q. 





\(a) K=2 ' 




- 


0.4 






- 


0.3 






- 


0.2 






n=o 


0.1 






-— -^oos 


0.0 


V/ 0.25 


(Xl>- 


__ai ' : 



0.5 



0.4 



0.3 



0.2 



0.1 - 



0.0 



V(b)K=2 


n=o/ \ 


IV 


1 1 ' 


I \ 0.5 

|\ X 0.25 


1 - 




~\ / \ 


, V^^r— - "-! i r i 


-V v 



50 60 10 20 30 40 50 60 70 80 90 100 

0.5 




50 60 20 40 60 80 100 120 140 160 180 200 



FIG. 3: Time evolution of the off-diagonal element P23 of the reduced density matrix of a Heisenberg two-spin system (J = —5), 
interacting with a Heisenberg-type environment He via an isotropic Heisenberg Hamiltonian Hse (A = —0.075 ) for the same 
geometric structures in the environment: (a,b) K — 2 and (c,d) K = N — 1. The number next to each curve is the corresponding 
value of fl. Adapted from Ref.— . 



quantum system is not sensitive to the presence of a 
Heisenberg-type environment He, but there is a decay of 
the amplitude of the long-living oscillations. The larger 
K (see Fig. C^b-e) (left)) or ft (see Fig.|3[a,c)), the faster 

the decay is. 

If |A| <C Qy/K and $7 is comparable with J, keeping 
K fixed and increasing $7 smoothly changes the initial 
decay from Gaussian (fast) to exponential (slow). The 
long-living oscillations are completely suppressed (see 
Fig. [Hb,d)). For large $7, the simulation data fits very 
well to 



\ P 2i(t)\= l -e- A «^\ 



(33) 



where Ak (O) is approximately linearly dependent on S7: 



A K (fi) » QA K , 



(34) 



and we find that Ai = 9.13 and An-i = 26.73. 

Physically, the observed behavior can be understood as 
follows 25 . If |A| rj fly/K, a bath spin is roughly equally 
affected by the motion of the other bath spins and the 
system spins. Therefore, each bath spin follows the orig- 
inal dynamics, as if there was no coupling between bath 



spins. This explains why the initial Gaussian decay is 
insensitive to the values of K or il. After the initial de- 
cay, the whole system is expected to reach a stationary 
state, but because of the presence of Heisenberg-type in- 
teractions between the bath spins, it leads to a decrease 
of the coherence between the singlet and triplet states, 
and therefore a new stationary state of the bath is estab- 
lished, suppressing the long-living oscillations. 

For larger K, the distance between two bath spins, de- 
fined as the minimum number of bonds connecting the 
two spins, becomes smaller. For instance, for K = 2, 
this distance is (N — 2) /2, and for K = N — 1, it is zero. 
For fixed $7 and larger K the fluctuations in the spin 
bath propagate faster, and therefore the evolution to the 
stationary state is faster. Furthermore, since the environ- 
ment in our model is a highly frustrated system, increas- 
ing the connectivity K will increase the energy resolution 
of the eigenstates, which makes the dynamics of the envi- 
ronment more complicated. For fixed K, increasing the 
coupling strength between the bath spins will speed up 
the dynamics of the bath, that is, the larger H, the faster 
will be the evolution to the stationary state. However 
the coupling strength within the environment should not 





20 40 60 80 100 120 140 160 180 200 



FIG. 4: Time evolution of the Loschmidt echo L (t) of a Heisenberg two-spin system (J = —5), interacting with a Heisenberg- 
type environment He via a Heisenberg (A = —0.075) Hamiltonian Hse- The values of Q,\fK are: (a) Q,\fK = O.lVAf — 1, 
(b) Q,\TK = 0.15y/W^T, (c) Q^K = 0.25ViV- 1, and (d) Q-/K = y/N - 1. The different lines in each pannel correspond 
to different if. Solid (black) line: K = 2; dashed (red) line: if = 4; dotted (green) line: K = 6, and dash-dotted (blue) line: 
if = AT - 1. Adapted from RefA 



be too large, because otherwise the energy resolution in 
the bath will be too small to lead the energy dissipation 
of the quantum system. 

In the opposite case |A| <C flVK and fl is compara- 
ble with J, Hse is a small perturbation relative to He 
and the coupling between the bath spins is the dominant 
factor in determining the dynamics of the bath spins. 
Therefore, by increasing K or Q, the bath spins will have 
less freedom to follow the dynamics induced by the cou- 
pling to the two system spins, the influence of the bath on 
the quantum system will decrease, and the (exponential) 
decay will become slower. 

Here we have compared Qy/K to |A| to distinguish 
different regimes. As a matter of fact, Qy/K does not 
completely characterize the decoherence process, but it 
can be used to characterize its time scale. Indeed, as 
shown in Fig. [H for different \[K and Q but the same 
value of ily/K, the time evolution of Lit) is very similar. 
Note that if fl\/K increases (compare Fig. Ht to Fig. EH), 
the differences between the Loschmidt echoes increase. 

According to the general picture of decoherence^, for 
an environment with nontrivial internal dynamics that is 



initially in a random superposition of all its eigenstates, 
we expect that the quantum system will evolve into a 
stable mixture of its eigenstates. In other words, the de- 
coherence will cause all the off-diagonal elements of the 
reduced density matrix to vanish with time. In the case 
of an isotropic Heisenberg coupling between the quan- 
tum system and the environment, H$ commutes with 
the Hamiltonian H, hence the energy of the quantum 
system is a conserved quantity. Therefore, the weight of 
the singlet \S) in the mixed state should be a constant 
(1/2), and the weights of the degenerate eigenstates |7b), 
|T_i) and \T{) are expected to become the same (1/6). 
As shown in Fig. [2jb-e) (right), our simulations confirm 
that this picture is correct in all respects. 



C. Decoherence With Energy Dissipation 

Now we consider the case that there is energy dissi- 
pation of the quantum system, i.e. [Hs,H] ^ 0. First, 
instead of considering a Heisenberg system-environment 
interaction Hse as in Fig. [5121 we now take a Heisenberg- 



10 



type Hse in Fig- [MSI the other interactions like Hs and 
He are the same as in Fig. [2E1 From a direct compar- 
ison of these results, it is clear that the roles of K and 
n are the same, no matter whether the energy of the 
quantum system is conserved or not. If |A| 3> fl\/K, 
in the presence of anisotropic interactions between the 
quantum system and the environment spins, the second 
step of the oscillations decay and finally disappear as 
K increases, even in the absence of interactions between 
the bath spins. This is because the anisotropic interac- 
tions break the rotational symmetry of the coupling be- 
tween the quantum system and the environment which 
is required for the long-living oscillations to persist. If 
|A| <C Q\/K and fl is comparable with J, |p 2 3 (t)\ can 
still be described by Eq. (|33|) . but now Ak (^2) is no 
longer a linear function of Q. This is because the energy 
dissipation will change the weight of each pointer state 
(eigenstate) in the final stable mixture, which makes the 
time evolution of |p23 (t)\ more complicated. 

More results with Isingzb interaction Hse are shown in 
Fig. [7] and Fig. [5] They give similar results as in the case 
of a Heisenberg-type interaction Hse- The environments 
in Fig. [5] are not a spin glass, but an isotropic Heisenberg 
antiferromagnetic square or triangle lattice, which are 
also frustrated systems. 



D. Summary 
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In conclusion, with a frustrated spin-bath environment 
that initially is in a random superposition of its basis 
states, a pure quantum state of the quantum spin sys- 
tem will evolve into a classical mixed state. If the in- 
teraction between the quantum system and environment 
is much smaller than the coupling between the spins in 
the quantum system, the pointer states are the eigen- 
states of the quantum system. Both of these observa- 
tions are in concert with the general picture of decoher- 
ence^. Furthermore, if the energy of the quantum system 
is conserved, the pointer states can still be the eigen- 
states of the quantum system, independent of the ratio 
of the system-bath coupling to the coupling within the 
system. For the anisotropic spin-bath, changing the in- 
ternal dynamics of the environment (geometric structure 
or exchange couplings) may change the decoherence of 
the quantum spin system from Gaussian to exponential 
decay. 



V. EVOLUTION TO THE GROUND STATE OF 
A TWO-SPIN SYSTEM 

To approach the ground state of a quantum system by 
coupling it to a large quantum bath is not a trivial prob- 
lem. The interactions between the two quantum systems 
will not only lead to the exchange of the energy, but also 
the coherence of the wave function. In this section, we 
show under which condition a two-spin system can ap- 



FIG. 5: Time evolution of the real part of the off-diagonal el- 
ement p23 (left panel) and the diagonal elements pn, . . . , P44 
(right panel) of the reduced density matrix of a Heisen- 
berg two-spin system (J = —5), coupled via an isotropic 
Heisenberg-type interaction Hse (A = —0.15 ) to a 
Heisenberg-type environment He (f2 = 0.15) with different 
connectivity: (a) K = 0; (b) K = 2; (c) K = 4; (d) K = 6; 
(e) K = N - 1. 



proach its ground or near-ground state. The affect of 
entanglement and interaction symmetry during this evo- 
lution will be discussed. 

The interaction within the two-spin system will be 
fixed as isotropic Heisenberg ferromagnet ( J > 0) or an- 
tiferromagnet (J < 0). In the case of ferromagnet, there 
are three degenerate ground states (triplet states) |T ), 
\T±) and |T_i), whereas in the case of antiferromagnet 
it is the singlet state \S). Therefore when the quantum 
system approaches the state which has the same energy 
as the ground state(s), it can be a single eigenstate if the 
ground state is non-degenerate, and it can also be an en- 
tangled superposition (quantum) or mixed state (classic) 
if the ground states are degenerate. 

In order to let the two-spin system approach the 
ground state, it is necessary to keep the environment 
at low temperature. Instead of a classical mixture, we 
will prepare the environment in a pure quantum state, 
in which the temperature is not well defined. So we 
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FIG. 6: Same as Fig. [3] except that Hse is Heisenberg-type and A = 0.15. Adapted from Ref.' 



simply initialize the environment in its ground or near- 
ground state to guarantee the one-direction energy flow. 
As in Rcf . 25 ' 26 , this initial ground state of the environ- 
ment leads to a sharp local density of states in the whole 
system, and therefore the decoherence of the quantum 
system is much weaker comparing to cases with random 
initial state of the environment. 

The interactions between the two-spin system and the 
spin bath will be set as Heisenberg-type or Ising-type. 
Even in both cases the interactions of different orienta- 
tions (x,y and z) are not all the same, but they are totally 
different: the Heisenberg-type Hse is still symmetrical 
because the exchange interactions have the same random 
amplitude in the three orientations, whereas the Ising- 
type Hse is antisymmetric since the exchange interac- 
tions are totally different between x(y) and z directions. 

The interactions within the environment are fixed as 
Heisenberg-type, which reduce the strength of the coher- 
ence between the quantum system and the environment, 
as we showed in the previous section and Refi^^r— . 



A. Symmetrical Coupling 

We first consider the case that an antiferromagnetic 
(J < 0) or ferromagnetic (J > 0) quantum system that 



interacts with the Heisenberg-type environment via a 
Heisenberg-type interaction. 

In Fig. [SJ we present simulation results for the two- 
spin correlation function (as a measure of the energy) 
for different values of the coupling strength (fi) in the 
environment. Clearly, in case (a), the relaxation in both 
cases, antiferromagnetic and ferromagnetic, is rather slow 
and confirming that there is relaxation to the ground 
state requires a prohibitively long simulation. For cases 
(b) - (d), the results are in concert with the intuitive 
picture of relaxation due to decoherence: The correlation 
shows the relaxation from the up-down initial state of the 
quantum system to the ground or near-ground state. As 
for the antiferromagnetic quantum system (J < 0), the 
two-spin correlation relaxes to a value of about 0.65 - 
0.70, which is much further away from the ground state 
value —3/4 than we would have expected on the basis of 
the results of the ferromagnetic quantum system. In the 
true ground state of the whole system, the value of the 
two-spin correlation in case (b) of J < is —0.7232, and 
hence significantly lower than the typical values, reached 
after relaxation. On the one hand, it is clear (and to be 
expected) that the coupling to the environment changes 
the ground state of the quantum system, but on the other 
hand, our numerical calculations show that this change 
is too little to explain the apparent difference from the 
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FIG. 7: Time evolution of the real part of the off-diagonal 
element p2z of the reduced density matrix of a Heisenberg 
two-spin system (J = —5), coupled via an Isingzb interaction 
Hse (A = 0.075 ) to a Heisenberg-type-ring environment He 
(N = 16) with different range of fi (0,0.01,0.05,0.1, 1,2,5). 
The value of the corresponding 57 is indicated by the number 
near each line. 



results obtained from the time-dependent solution. In 
case (e), the characteristic strength of the interactions 
between the spins in the environment is of the same order 
as the exchange coupling in the quantum system (57 « J), 
a regime in which there clearly is significant transfer of 
energy back-and-forth, between the quantum system and 
the environment. 

In Fig. [TU1 we show the diagonal elements of the re- 
duced density matrix for case (b). After reaching the 
steady state, the nondiagonal elements exhibit minimum 
fluctuations about zero and are therefore not shown. 
From Fig. [TU1 it is then clear that the quantum sys- 
tem relaxes to the singlet state for the antiferromagnetic 
(J < 0) system, and to a mixed state for the ferromag- 
netic system (J > 0), as expected on intuitive grounds. 

An important observation is that our data convincingly 
shows that it is not necessary to have a macroscopically 
large environment for decoherence to cause relaxation to 
the ground state: A spin-glass with N = 14 spins seems 
to be more than enough to mimic such an environment 
for a two-spin system. This observation is essential for 
numerical simulations of relatively small systems to yield 
the correct qualitative behavior. 

Qualitative arguments for the hig h efficiency of the 
spin-glass bath were given in Ref. J2J- Since spin-glasses 
possess a huge amount of states that have an energy close 
to the ground state energy but have wave functions that 
are very different from the ground state, the orthogo- 
nality catastrophe, blocking the quantum interference in 
the quantum system^ is very strongly pronounced in 
this case. 

From the data for (b) - (d), shown in Fig. [9l we con- 
clude that the time required to let the quantum system 



relax to a state that is close to the ground state depends 
on the energy scale (57) of the random interactions be- 
tween the spins in the environment. As it is difficult to 
define the point in time at which the quantum system 
has reached its stationary state, we have not made an 
attempt to characterize the dependence of the relaxation 
time on 57. 



B. Antisymmetric Coupling 

In our simulation, the initial state of the quantum sys- 
tem is |tl) and this state has total magnetization M = 0. 
For an Ising-type interaction Hse of the envionment with 
a Heisenberg system Hs, the magnetization M of the 
quantum system commutes with the Hamiltonian of the 
whole system. Therefore, the magnetization of the quan- 
tum system is conserved during the time evolution, and 
the quantum system will always stay in the subspace with 
M = 0. In this subspace, the ground state for the an- 
tiferromagnetic quantum system is the singlet state \S) 
while for the ferromagnetic quantum system the ground 
state (in the M = subspace) is the entangled state \T ). 
Thus, for the Ising-type interaction Hse, starting from 
the initial state |f4-)i the quantum system should relax to 
an entangled state, for both a ferro- or antiferromagnetic 
quantum system, that is, at any time t, the state of the 
whole system can be written as 



|*(t)) = |S}|$ s (i)) + |T }|$ To (i)), (35) 

where |$s) and |$t d ) denote the states of the environ- 
ment. 

Let us denote by {|<fo)} a complete set of states of the 
environment. Within the subspace spanned by the states 
{|,S)|^>i), |Xo)|0i)}j the Hamiltonian can be written as 

H = E s \S)(S\+E t \T )(T \+He 

N 



1 



A*h 



-£(Ag-A^)(|S)(T | + |T )(S|)/;, 

(36) 

where we used (S\Sf\S) = (T \Sf\T ) = (S\S%\S) = 
(To\SI\T ) = 0, (T \Sf\S) = 1/2, and (T \S$\S) = -1/2. 
Introducing a pseudo-spin a = (a x , a v ,a z ) such that 
the eigenvalues +1 and —1 of a z correspond to the states 
\S) and \Tq), respectively, Eq. (151)1) can be written as 

E s - E T Es + Et „ 

ti = 1 a + He 



1 N 



(37) 



j=i 



showing that in the case of Ising-type Hse, the quantum 
system with two spins is equivalent to the model Eq. (j37p 
with one spin. 
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FIG. 8: Time evolution of the real part of the off-diagonal element p2z and concurrence of the reduced density matrix of 
a Heisenberg two-spin system (J = —1), coupled via an Ising interaction ± Hse (A = 0.075 ) to an isotropic Heisenberg 
environment He (fi = —0.15, N — 16) with connectivity K = 4 (square lattice) and K = 6 (triangle lattice). 



From Eq. (|3T|) . it follows immediately that the Hamil- 
tonian is invariant under the transformation {J,a z } — > 
{ — J, — a z }. Indeed, the first, constant term in Eq. ((37)) is 
irrelevant and we can change the sign of the second term 
by rotating the speudo-spin by 180 degrees about the x- 
axis. Therefore, if the initial state is also invariant under 
this transformation, the time-dependent physical proper- 
ties will not depend on the choice of the sign of J, hence 
the ferro- and antiferromagnetic system will behave in 
exactly the same manner. In our case, the initial state 
can be written as (\S) + |Tb))|e/>o)/\/2, which is trivially 
invariant under the transformation a z — > —o z . 



Therefore for Ising-type Use (A 2 - 



Av) 



0), an 



initial state that is invariant for the transformation \S) O 
|T )), (^(t)\A\^f(t)) does not depend on the sign of J, for 
any observable A of the quantum system that is invariant 
for this transformation. Under these conditions, it is easy 
to prove that 



(*(t)|S 1 -S 2 |*(t))f + (*(t)|Si 



S 2 Mt))A = -\, (38) 



where the subscript F and A refer to the ferro- and an- 
tiferromagnetic quantum system, respectively. 

Likewise, for the concurrence we find Cf (t) = Ca (i) 
and similar symmetry relations hold for the other quan- 
tities of interest. Of course, this symmetry is reflected in 
our numerical data also, hence we can limit ourselves to 
present data for the antiferromagnetic quantum system 
with Ising-type Hse- 

In Fig.rTl"land ll2l we present simulation results for the 
two-spin correlation function and concurrence for differ- 



ent values of the parameter A and f2. It is clear that 
for a certain range of the interaction strength, the quan- 
tum system relaxes to a state that is very close to the 
ground state, see Fig.fTTtc) and[T2Ic). That is, the pres- 
ence of a conserved quantity (the magnetization of the 
quantum system) acts as a catalyzer for relaxing to the 
ground state. Intuitively, we would expect that the pres- 
ence of a conserved quantity hinders the relaxation and 
indeed, the relaxation in Fig. [TT] is much slower than in 
Fig. [9] Notwithstanding this, in the presence of a con- 
served quantity, the quantum system relaxes to a state 
that is much closer to the true ground state than the 
one it would relax to in the absence of this conserved 
quantity. 



One should notice that only in a small range of param- 
eters A and $7 the two-spin quantum system can evolve 
into a near-ground state. If the interactions within the 
bath is too small, the range of the energy spectrum of 
the bath limits the energy dissipation of the system, and 
the system spins mainly follow its own dynamics just 
like there is no bath. On the contrary, if the interactions 
within the bath are too strong, a small change of the con- 
figuration of the bath spins will lead to a large change of 
its energy, therefore the direction of the energy flow will 
oscillate with time and the quantum system cannot arrive 
to a state with approximately stable energy. 
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FIG. 9: Time evolution of the correlation (*(f)|Si-S 2 |*(t)) of 
the antiferromagnetic (top panel) and ferromagnetic (bottom 
panel) two-spin quantum system with Heisenberg-type inter- 
action Hse and environment He- The model parameters are 
A = 0.15 and a: Q = 0.075; b: Q = 0.15; c: Q = 0.20; d: 
fl = 0.30; e: fl = 1. The number of spins in the environment 
is N = 14. 



C. Summary 

In general, it turns out that the relaxation to the 
ground state is a more complicated process than one 
would naively expect, depending essentially on the ratio 
between parameters of the interaction and environment 
Hamiltonians. Two general conclusions are: 

(1) the quantum system more easily evolves into its 
ground state when the latter is more degenerate (degen- 
erate triplet states compared to the singlet) or less en- 
tangled (e.g. up-down state compared to the singlet); 

(2) constraints on the system such as existence of addi- 
tional integrals of motion can make the evolution to the 
ground state more efficient. 

An explanation of the first statement is that if the 
ground state is more degenerate, its affective dimen- 
sion in the Hilbert space is larger and therefore the 
wave function will have more possibilities to evolve into 
this subspace, especially if the environment is frustrated. 
The latter statement looks a bit counterintuitive since it 
means that it may happen that a more regular system 



FIG. 10: Time evolution of the diagonal matrix elements 
of the reduced density matrix of the antiferromagnetic (top 
panel) and ferromagnetic (bottom panel) two-spin quantum 
system for A = 0.15 and fi = 0.15 (case (b) of Fig. [9] except 
that the number of spins in the environment is N — 16) 



exhibits stronger relaxation than a chaotic one. The rea- 
son that it may happen is that introducing an additional 
integral of motion, such as the total magnetization, limits 
the dimensionality of the available Hilbert space for the 
quantum system. The larger the dimensionality of the 
available Hilbert space, the more complicated the deco- 
herence process is due to the appearance of the whole 
hierarchy of decoherence. A manifestation of this phe- 
nomenon has been observed earlier—: Under certain con- 
ditions, the same quantum system as studied here (four 
by four reduced density matrix) displays "quantum oscil- 
lations without quantum coherence" whereas for a single 
spin in a magnetic field (two by two reduced density ma- 
trix) decoherence can, relatively easily, suppress the Rabi 
oscillations completely 



VI. DECOHERENCE OF A MANY-SPIN 

SYSTEM 

In the previous sections, we focused on the quantum 
system with only two spins. Starting from this section 
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FIG. 11: Time evolution of the correlation (*(t)|Si •S 2 |*(t)> 
of the antiferromagnetic quantum system with Ising-type 
Hse and Heisenberg-type He- The model parameters are 
(a-d) A = 0.075, (e) A = 0.0375, (f) A = 0.15, and (a) 
fi = 0.075, (b,e) fi = 0.15, (c,f) Q = 0.30, (d) Q = 1. The 
number of spins in the environment is N — 16. 



we will consider more complicated systems which contain 
more spins, e.g. four to eight spins. 

In general when a quantum system interacts with a 
quantum environment, there will be energy dissipation 
and entanglement of their wave functions. This entangle- 
ment does not necessary lead to a classical mixed state 
of the quantum system, especially if the environment has 
finite size. But introducing certain properties, a finite 
quantum environment can also drive the quantum sys- 
tem to an exact classical state. In the previous sections, 
we showed that a frustrated environment can enhance 
the decoherence of a two-spin system, and as we will see 
in what follows, this is also the case for the many-spin 
system. 



A. Origin of the Microcanonical Distribution 



sity matrix of the system is calculated every r = 7r/10. 
In both cases, the diagonal terms of the reduced density 
matrix of the Heisenberg-ring approach a stable value af- 
ter an initial decay or increase, and the off-diagonal terms 
are all zero (a — > 0), which means that in both cases, the 
four spin quantum system approaches a fully decoherent 
state. 

Another significant result is that the degenerated en- 
ergy eigenstates have the same weigth distribution in the 
fully decoherent state, i.e. 7 — > 0, indicating that the fi- 
nal system is a microcanonical state in each eigenenergy 
subspace. The diagonal elements of the reduced density 
matrix corresponding to degenerate eigenstates with zero 
weight in the initial state are zero and remain so during 
the time revolution, see P5-11 in Fig. 1131 This is sim- 
ply due to the conservation of energy of the system, and 
due to the fact that the time evolution operator e~ tHt 
prevents changing the coefficients corresponding to these 
eigenstates. 

In Fig. Q3] we plot the time evolution of a (t) with 
different random realizations for the initial state of the 
environment or for the model parameters f2f and A"-. 
The difference between the curves is very small, indicat- 
ing that in our model, a particular randomness of the 
coupling parameters or initial state is not relevant to the 
general properties of the simulation results. 

In Fig. QJ3 we show the time evolution of a (t) for dif- 
ferent coupling strengths (O) in the environment. In gen- 
eral, increasing the coupling strength within the environ- 
ment will increase the effective energy range of the bath, 
which leads the decoherence more completely. But as 
we have shown in the case of the two-spin system, the 
coupling strength should not be too large, otherwise the 
energy resolution of the bath will be too small to lead 
the full decoherence of the system. 

In Fig. [Tni we show the time evolution of lg[cr (<)] for 
the same systems but with different topological struc- 
tures (connectivity K) in the environment. It is clear 
that as soon as there is frustrated interaction within the 
environment (K > 0), no matter what kind of topological 
structure it is, the decoherence of the quantum system is 
quite similar. 



Let's first consider the case that there is no energy dis- 
sipation or the energy dissipation is very small. In this 
case the decoherence of the quantum system fully origi- 
nates from the phase correlation of the environment. In 
Fig. [T21 a Heisenberg-ring with four spins is coupled to a 
frustrated spin glass consisting of 14 spins, and the inter- 
action between the quantum system and the environment 
is isotropic or anisotropic but very small (comparing to 
the coupling within the system, i.e. A = 0.15 <C \J\ = 5). 
In the former case the energy of the quantum system is 
conserved, and in the latter case the energy dissipation is 
very small so that it can be ignored. The system has four 
distinct eigenvalues (E\ = —2, E 2 -i = — 1, £5-11 = 0, 
and -E12-16 = 1) and sixteen different eigenstates. Dur- 
ing the time-integration of the TDSE, the reduced den- 



B. Summary 

If there is no energy dissipation, or the energy dissi- 
pation is so small that it can be ignored, then the en- 
tanglement between the quantum system and the envi- 
ronment occurs only in the subspace of the (degenerate) 
eigenstates which have nonzero weigth distribution in the 
initial state of the quantum system. That is, the possible 
pointer states in the mixed states are determined by the 
initial state of the quantum system itself. If a particular 
environment can lead to the decoherence of the quan- 
tum system without energy dissipation, then turning on 
the energy dissipation will still lead to decoherence, and 
even more completely. In fact, energy dissipation is not 
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FIG. 12: Time evolution of the concurrence of the antiferromagnetic quantum system with Ising-type interaction Hse and 
Heisenberg-type environment He- The model parameters are A = 0.075 and (a) Q = 0.075, (b) fi = 0.15, (c) £2 = 0.30, (d) 
= 1. The number of spins in the environment is N — 16. 



related to the question whether a quantum system can 
evolve into a classical mixed state or not. The main dif- 
ference between decoherence with or without energy dis- 
sipation are the number of possible pointer states in the 
mixed state. That is, an additional integral of motion of 
the system will limit the number of pointer states, and 
therefore a full decoherence state is a mixture with micro- 
canonical distribution in each eigenenergy subspace un- 
der the extra conservation law. On the contrary, if there 
is enough energy dissipation between the two systems, 
then all eigenstates of the quantum system are possible 
pointer states. And more importantly, as we will show in 
the next section, the mixed state of the quantum system 
follows the canonical distribution^. 



VII. THERMALIZATION OF A MANY-SPIN 
SYSTEM 

In the previous section, we have shown that turning 
on the interaction between the many-spin system and 
the environment, leads to a reduction of the coherence 
in the quantum system. The coupling with the environ- 
ment causes the initial pure state of the quantum system 



to evolve into a mixed state, obtained by tracing out all 
the degrees of freedom of the environment. The pointer 
states in the mixed state are determined by the initial 
state of the quantum system if there is an additional 
integral of motion, e.g.. the conserved energy or magne- 
tization. This leads to a microcanonical ensemble under 
a certain conservation law. On the other hand, if there 
is enough energy dissipation without any additional in- 
tegral of motion, we expect that the mixed state is a 
canonical ensemble^,. 

Earlier demonstrations of the fact that the system can 
be in the canonical ensemble state are based on Ergodic 
average d or canonical typicality 47-53 . The Ergodic 

averages consider the dynamics of a closed quantum sys- 
tem, and prove that in certain quantum systems, the 
expectation values of the dynamical variables of the sys- 
tem approach their values for the subsystem that is in 
the thermal equilibrium state. This is similar to the as- 
sumption of classical statistical physics, that is, during 
large enough time, the trajectory of the many particle 
system in the phase space will pass all possible points, 
and therefore the average of these points in the phase 
space will follow a certain distribution. On the contrary, 
the canonical typicality does not consider the dynamics of 
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FIG. 13: Time evolution of the diagonal terms (pi) and sum of the absolute values of the off-diagonal terms (a) in the reduced 
density matrix of a Heisenberg-ring Hs (J = —5, ns — 4, initial state \UD)s) coupled to a spin glass environment He (0 = 0.15, 
n = 16, initial state \RANDOM)e) via (a) Heisenberg interaction Hse (A = 0.075) or (b) Heisenberg-type interaction Hse 
(A = 0.15). Full decoherence is observed in both cases, and the system S relaxes to a state with equal weights within each 
energy subspace, that is, 7 — > 0, a microcanonical ensemble. 
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FIG. 14: The system represented here is the same as the one 
in Fig. I13f a). except that the randomness in the coupling 
constants of Hse or He, or the randomness in the initial 
state of the environment, are different in the curves A — E. 
It is clear that the time evolution of a(t) is not sensitive to 
the different random values of the coupling constants or the 
initial state of the environment if they follow the same type 
of random distribution. 



the system but assumes that if the whole system is in the 
microcanonical ensemble then its subsystems are in the 
canonical ensemble. The statement of the canonical typ- 
icality is quite general, but it is clear that the dynamical 
procedure is missing. Moreover, the coupling between the 
quantum system and the environment is assumed to be 
so small that it can be neglected in the theory. Therefore 
it is important to show how the canonical distribution of 
statistical mechanics relates to the dynamical evolution 
of the quantum system. 



A. Origin of the Canonical Ensemble 

As a frustrated environment is very effective for creat- 
ing full decoherence (er — > 0) in a quantum spin system, 
and full decoherence is a necessary condition for the state 
of the system to converge to its canonical distribution, 
we have chosen spin glass environments, which have no 
obvious symmetries 26 . 

First, we consider a system (Heisenberg-ring Hg) in- 
teracting (Heisenberg-type Hse) with an environment 
(spin glass He)- The system has four distinct eigenvalues 
(Ei = -2, £ 2 -4 = -1, £5-11 = 0, and E l2 -i6 = 1) and 
sixteen different eigenstates. The environment has 2 18 
eigenstates. During the time-integration of the TDSE, 
the reduced density matrix of the system is calculated ev- 
ery r = 7r/10 as in the previous section As we described 
earlier, the values of the diagonal elements pa yield an 
estimate for the effective inverse temperature b(t), the 
error S(t) for this estimate and the measure a(t) for the 
deviation from a non-diagonal matrix. The energy of the 
system is obtained as Es(t) = Trsp(t)Hs- 

From the simulation results, shown in Fig. II Tl it is 
clear that for t > 50r, each diagonal element pu of the 
reduced density matrix converges to one out of four sta- 
tionary values, corresponding to the four non-degenerate 
energy levels of the system. This convergence is a two- 
step process. First the system looses all coherence, as 
indicated by the vanishing of a (t) for t > 50r. The time 
dependence of a (t) fits very well to an exponential law 



a(t) 



Ae 



-t/T 2 



(39) 



with o-oo = 0.00128, A = 0.602 and T 2 = 8.01r. In 
the small panels of Fig. [TT1 the red dots are the simula- 
tion data and the black curves are the fitting function. 
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FIG. 15: Same as Fig. I13f a) except that the range of the coupling strength (fi) in the environment is different. 




and 



FIG. 16: Same as Fig.[l3ja) except that the topological struc- 
ture (connectivity K) of the environment is different. 



Likewise, the vanishing of S(t) on the same time-scale 
(T2 = 7.32t), indicates that the density matrix of the 
system converges to the canonical distribution with the 
same speed of decoherence. 

The effective temperature b(t) and the energy of the 
system E$ (t) also fit very well to the exponential laws 



E (t) =E 00 +Ce 



-t/T! 



(41) 



b(t)=(3 + Be-*/ Tl 



(40) 



with /3 = 0.0962, B = -0.900, and T x = 13.3r and 
Eoo = -0.0745, C = -0.952. The estimated values for 
T\ and Ti change very little if we choose different random 
realizations for the initial state of the environment or for 
the model parameters fi"- and A"- (see Fig. [TS]). If we 
change their range, T\ and T2 also change, as naively 
expected. 

In order to verify the role of the dynamics within the 
bath to the thermalization of the quantum system, we 
plot the time evolution of 6 (t) for different coupling 
strength (f2) in Fig. Q1JJ Similar as the dependence of a 
on fl in Fig. [15j increasing the coupling strength (within 
a certain range) will increase the effective energy range 
of the bath, which leads to a more complete decoherence 
and thermalization. 

The simulation results of a similar system with one 
extra spin in the quantum system (n s = 5) and one spin 
less in the environment (n — 17) are shown in Fig. 1211 
and are very similar to the ones shown in Fig. 1171 These 
simulations demonstrate that the system first looses all 
coherence and then, on a longer time-scale, relaxes to its 
thermal equilibrium state with a finite temperature. In 
terms of the theory of magnetic resonance^, T\ and T-i 
are the times of dissipation and dephasing, respectively. 
In the case of very small He, one should expect, instead 
of an exponential decay of a and E, a Gaussian decay, as 
observed in Ref<22r— . 



19 




5' 6'-?' 8 '-9 no ni 12 " u 15 1B 

J I I I I I I I I I I I I I I I I L_ 



80 100 120 140 160 

t/t 

0.6- 



180 200 




100 



FIG. 17: Simulation results for the diagonal elements p* = pu(t) of the density matrix of S, the energy Es = Es{t), the effective 
inverse temperature b = b(t) and its variance 8 = S(t), and a = a(t) which is a measure for the decoherence in S, as obtained 
by solving the TDSE for the whole system with a Heisenberg-ring Hs (J = — 1, «s = 4), a Heisenberg-type interaction Hse 
(A = 0.3), a spin glass environment He (CI — 1, n — 18), and r = 7r/10. The initial state of the whole system is a product 
state \UD) s <g> | RANDOM) E . The red dots in the small panels represent the simulation data, and the black curves are fitting 
curves (see text). 



It is necessary to extend the types of interaction 
(Hamiltonian) to verify the generality of the above re- 
sults. The time evolution of 8 (t) for the same systems 
but with different topological structures (connectivity K) 
in the environment are shown in Fig. 1201 Increasing the 
connectivity in the environment will increase the energy 
resolution and make the dynamics in the environment 



more complicated, which enhances the approach to the 
canonical distribution of the system. This is quite dif- 
ferent from the effect of the topological structure on the 
decoherence shown in Fig. [TH where as soon as there is 
frustrated interaction within the system (K > 0), the 
decoherence of the quantum system is quite similar for 
different topological structures. 



20 



eo. 



1.2 
1.0 
0.8 
0.6 
0.4 
0.2 



\ 



-A 
B 
C 

D 

E 



10 20 30 40 50 60 70 80 90 100 
t/T. 



FIG. 18: Same as Fig. 1171 except that the randomness in the 
coupling constants of Hse or He, or the randomness in the 
initial state of the environment, are different in the curves 
A — E. It is clear that the time evolution of the effective tem- 
perature b(t) is not sensitive to the different random values of 
the coupling constants or the initial state of the environment 
if they follow the same type of random distribution. 





FIG. 20: Same as Fig. 1 171 except that the topological structure 
(connectivity K) of the environment is different. 
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FIG. 19: Same as Fig. [17] except that the range of the coupling 
strength (fi) is different in the environment of each curve. 



FIG. 21: Same as Fig. 1171 except that ns = 5 and n = 17. 



More results for the system with different symmetries 
and connectivities but with the same type of environ- 
ments (He) and the same type of interactions (Hse) are 
shown in Fig. [22j The systems used are a XY-ring, a 
Heisenberg-ring, an Ising-ring, a Heisenberg-triangular- 
lattice, and a spin glass. From the results represented 
in Fig. [2U it is clear that independent of the internal 
symmetries and the connectivity of the system, and in- 
dependent of the initial state of the whole system, all 
systems relax to a state with full decoherence, except 
case f. The main difference between case f and the other 
cases is the initial state of the environment. Since the 
environment in our model is a highly frustrated system, 
the LDOS of the whole system covers the whole spec- 
trum of the energy space, no matter the initial state of 
the environmental spins is all spins up or all spins down, 
or random spins up and down, or a random superposi- 



tion of all the states in the spin up-down basis. But if 
the environment is prepared in the ground state or near 
ground state, then the LDOS of the whole system be- 
comes more sharply peaked, see b and f in Fig. [23l Up 
to a trivial normalization factor, the LDOS curve for case 
b is indistinguishable from the density of states (data not 
shown) calculated from the solution of the TDSE using 
the technique described in Rcf. 88 . This suggests that if 
the environment starts from the random superposition of 
all its energy eigenstates, all states of the whole system 
may participate in the decoherence/relaxation process. 
In contrast, the LDOS curve for case f has a very small 
overlap with the density of states. Therefore, starting 
with an environment in the ground state, only a rela- 
tively small number of states participates in the decoher- 
ence process, as confirmed by the results for a(t) shown 
in Fig. |52f. 

One should also notice that in case b, a vanishes expo- 
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FIG. 22: Simulation results for the energy Es = Es{t), the effective inverse temperature b = b(t), its variance S = S(t), and 
the deviation from a diagonal matrix a = a(t) as obtained by the solution of the TDSE for a variety of different systems 
S coupled to a spin glass environment He via a Heisenberg-type interaction Hse- The systems used are a: XY-ring, b 
and f: Heisenberg-ring, c: Ising-ring, d: Heisenberg-triangular-lattice, and e: spin glass. The initial states of the whole 
system are a: \GROUND) s (S \ RANDOM) E , b: \UD) S <g> \RANDOM} B , c: \UU)s <8> \RR)e, d: \UU)s <8> \RANDOM) E , 
e: \GROUND) s ® |t/Z3)_E, and f: |t/D) s <g> \GROUND)e- The numbers of spins in the system are ns — 8 for cases a-c and 
ns = 6 for cases d-f. The number of spins in the environment is n = 16 for all cases. The model parameters are J = — 1, 
A = 0.3 and Q — 1, except for case e in which A = 1. 



nentially with time, whereas in the other cases (a,c,d,e), 
a initially increases and then vanishes exponentially with 
time, due to the entanglement between the system and 
the environment. This observation is in concert with our 
earlier worh^S 5 .. Furthermore, in all cases except f, the 
system always relaxes to a canonical distribution (<5 — > 0) 
as soon as it is in the state with full decoherence (<j — >• 0), 
indicating that the time of decoherence (T2) and the time 
of thermalization (Ti) is almost the same. In agreement 
with the results depicted in Fig.[17l the decoherence time 
T2 is shorter than the typical time scale T\ on which the 



system and the environment exchange energy. Note that 
in contrast to the cases considered in the theory of nuclear 
magnetic resonance, in most of our simulations, Hs, He 
and Hse are comparable so the standard perturbation 
derivation of a and E does not work. 

The negative temperature (b < 0) is also observed in c 
and d. In fact, as the temperature T is defined as 

- - — (42) 

T dE 1 [ ' 

where S and E are the entropy and energy. In a quantum 
spin system, the entropy may decrease when the energy 
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FIG. 23: Simulation results for the local density of states 
as a function of the energy. Solid line: Case corresponding 
to Fig. [23j. The initial state is \UD) S ® \RANDOM) E \ 
Dashed line: Case corresponding to Fig. I22f . The initial state 
is \UD) S ® \GROUND) E . Adapted from RefA 



increases. For example, the states with all spins up and 
all spins down have the same entropy, but can have to- 
tally different energy Suppose there is no interaction be- 
tween N ferromagnetic spins (J > 0), and there is a uni- 
form magnetic filed applied on the Z-\- direction, then the 
state with all spins up is the ground state, and with all 
spins down is the eigenstate with highest energy. Chang- 
ing the magnetization of the states from M max = N/2 to 
Mmin = —N/2 will change the sign of the temperature 
T(M) at the point M = 0, that is, T > when M > 
and T < when M < 0. 



B. Summary 

We have shown that if we have a system that inter- 
acts with an environment and the whole system forms a 
closed quantum system that evolves in time according to 
the TDSE, then a frustrated environment with a random 
distribution in the energy basis will lead to the full de- 
coherence of the system. Furthermore, if the system and 
environment can exchange energy, the range of energies 
of the environment is large compared to the range of en- 
ergies of the system, then the mixed state of the system 
is a canonical distribution. 



VIII. CONCLUSION AND DISCUSSION 

The results presented here have been obtained from an 
ab initio numerical solution of the TDSE in the absence 
of, for instance, dissipative mechanisms, and demonstrate 
that the existence of the microcanonical distribution (in 
each eigenenergy subspace) and the canonical ensemble 
is a direct consequence of quantum dynamics. 



We emphasize that our conclusion does not rely on 
time averaging of observables, in concert with the fact 
that real measurements of thermodynamic properties 
yield instantaneous, not time-averaged, values. Further- 
more and perhaps a little counter intuitive, our results 
show that relatively small environments (~ 20 spins) are 
sufficient to drive the system (2 — 8 spins) to thermal 
equilibrium and that there is no need to assume that 
the interaction between the system and environment is 
weak, as is usually done in kinetic theory. Note that even 
if most cases shown in this review clearly indicate a full 
decoherence of the quantum system, it does not mean 
that the condition of full decoherence is a nontrivial re- 
quirement. 

To conclude, we find that: 

(1) Frustration of the interactions, either within the 
environment or between the system and the environment, 
enhances the decoherence of the quantum system. 

(2) The quantum system more easily evolves into its 
ground state when the latter is more degenerate or less 
entangled, or has certain additional integrals of motion. 

(3) The distribution of the state of a quantum system 
is the microcanonical or canonical ensemble only if the 
system is in a fully decoherent state. 

(4) The restriction of a fully decoherent state to be a 
microcanonical ensemble per eigenenergy subspace is the 
presence of an additional integral of motion except a con- 
served energy. For example, a conserved magnetization 
of the quantum system prevents parts of the degenerate 
eigenstates to be the pointer states. 

(5) The restriction of a microcanonical ensemble to 
be a canonical ensemble is the presence of an additional 
integral of motion, so that the energy of each subspace is 
conserved. 

(6) The distributions in quantum statistical mechanics, 
such as the microcanonical and canonical distributions, 
are the direct consequence of quantum dynamics. 

Finally we want to discuss the second law of thermo- 
dynamics in quantum systems. 

The second law of thermodynamics states that the en- 
tropy of an isolated system which is not in equilibrium 
tends to increase over time, approaching a maximum 
value at equilibrium. In quantum mechanics, the state 
of a closed quantum system is always a pure state and 
therefore its entropy is a constant (zero). It is thus clear 
that the second law of thermodynamics is not valid in a 
closed quantum system. If a quantum system starts to 
interact with an environment, its entropy increases from 
zero but may not reach a maximum value at equilibrium. 
The dynamics of the whole system could be periodical 
and therefore the time evolution of the states could be 
reversible. 

Our numerical results show that if the quantum sys- 
tem becomes a classical mixed state, then the time evolu- 
tion becomes irreversible and the entropy becomes stable 
when it reaches the maximum value. In fact, the state 
with the microcanonical distribution in each eigenenergy 
subspace or in the canonical ensemble has the maximum 
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entropy within all possible states that the system could 
be. For a quantum system with a certain energy and 
a fixed number of particles it is such that, if there is a 
conservation law to restrict some eigenstates in the re- 
duced density matrix, then the state with maximum en- 
tropy corresponds to the one with all the accessible de- 
generate states having the same weigth distribution, i.e. 
a microcanonical distribution per eigenenergy subspace. 
On the other hand, if there is no such restriction on the 
eigenstates, then the canonical ensemble is the state with 
maximum entropy, as proved by Jaynes in Ref»22. 

We may conclude that the validity of the second law 
of thermodynamics in quantum mechanics is related to 
the decoherence process of the quantum system. If a 
quantum system becomes classical under the influence 
of the environment, then its entropy will increase until it 



reaches a maximum value of all possible mixed states, i.e. 
the microcanonical distribution per eigenenergy subspace 
or the canonical ensemble. If a quantum system cannot 
evolve into a stable mixed state, then its entropy will not 
always follow the second law of thermodynamics. 
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